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Abstract. We prove that if p > 1 and < r < p then the sequence ( m ^ r )^+^, 
to = 0,1,2,..., is positive definite, more precisely, is the moment sequence of a prob- 
ability measure n(p,r) with compact support contained in [0, +oo). This family of 
measures encompasses the multiplicative free powers of the Marchenko-Pastur distri- 
bution as well as the Wigner's semicircle distribution centered at x = 2. We show 
that if p > 1 is a rational number, < r < p, then /x(p, r) is absolutely continuous 
and its density W P:r (x) can be expressed in terms of the Meijer and the generalized 
hypcrgcometric functions. In some cases, including the multiplicative free square and 
the multiplicative free square root of the Marchenko-Pastur measure, W Ptr (x) turns 
out to be an elementary function. 
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Introduction 

For p, r e K we define the Raney numbers (or two-parameter Fuss-Catalan numbers) 



m—l 

r 



by 

(1) A m (p,r) := -^r JJ(mp + r - i), 

m - i= i 

A (p, r) :— 1. For m — 0, 1, 2, ... we can also write 

'mp + r\ r 
m J mp + r ' 

(provided mp + r ^ 0) where the generalized binomial is defined by 



(2) A m (p,r) = 



a \ a(a — 1) . . . (a — m + 1) 



mj mi 

Let B p (z) denote the generating function of the sequence {A m (p, 1)}^ =0 , the Fuss num- 
bers of order p: 

oo 

(3) B p (z):=J2A m (pA> m , 

m=0 

convergent in some neighborhood of 0. For example 

(4) B 2 (z) = ' 



l + ^l^Iz" 
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Lambert showed that 

oo 

(5) B p (zy = Y,A m (P,r)z m , 

m=0 

see [8] . These generating functions also satisfy 

(6) B p (z) = 1 + zB p (zf, 
which reflects the identity A m (p,p) = A m+ i(p, 1), and 

(7) B p (z)=B p . r (zB p (z) r ). 

It was shown in [13] that if p > 1 and < r < p then the sequence {A m (p, r)}^ =0 
is positive definite, i.e. is the moment sequence of a probability measure fi(p, r) on 
M. Moreover, /i(p, r) has compact support (and therefore is unique) contained in the 
positive half-line [0, oo) (for example fi(p, 0) = 5o). The proof involved methods from 
the free probability theory (see (231 C3 E])- I n particular, for p > 1 

(8) A i(P,l)=M2,lf (p " 1) , 

where Kl denotes the multiplicative free power, and /z(2, 1) is known as the Marchenko- 
Pastur (called also the free Poisson) distribution. It is given by 



(9) M (2,1) = ^W — dx on [0,4], 

Z7T V X 

and plays an important role in the theory of random matrices, see (24"l [9j [101 El HI H] • It 
was proved in [1] that the measure /z(2, l) Hn = /i(n + l, 1) is the limit of the distribution 
of squared singular values of the power G n of a random matrix G, when the size of the 
matrix G goes to infinity. 

In this paper we are going to prove positive definiteness of {A m (p, r)}^ =0 using more 
classical methods. Namely, we show that ifp>l,0<r<p and if p is a rational number 
then fi(p, r) is absolutely continuous and can be represented as Mellin convolution of 
modified beta measures. Next we provide a formula for the density W p>r (x) of n(p,r) 
in terms of the Meijer function and consequently, of the generalized hypergeometric 
functions (cf. (251 HE], where p was assumed to be an integer). This allows us to 
draw graphs of these densities and, in some particular cases, to express W PtT (x) as an 
elementary function. It is worth to point out that for r = 1 an alternative description of 
the densities W p i(x) has been recently given by Haagerup and Moller, see Corollary 3 
in [IT]. 

Finally let us also mention that the measures /i(p, r) satisfy a peculiar relation: 

(10) fi(p,r)> fi(p + s,s) = [i(p + s,r + s), 

for p > 1, < r < p and s > 0, see [13] . involving monotonic convolution ">", an 
associative, noncommutative operation on probability measures on R, introduced by 
Muraki [H]. 

1. Preliminaries 

For probability measures fix, fi2 on the positive half-line [0, oo) the Mellin convolution 
is defined by 

poo poo 

(11) (nx o /j, 2 ) (A) := / / l A (xy)dfi 1 (x)dfi 2 (y) 

Jo Jo 
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for every Borel set A C [0, oo). This is the distribution of product X\ ■ X 2 of two 
independent nonnegative random variables with Xi ~ /Xj. In particular, if c > then 
fj, o S c is the dilation of /i: 

(p5 c ) (A)=B cfJ i{A) :=fi(-A). 



Note that if /i has density /(x) then D c (/x) has density f(x/c)/c. 
If both the measures /ii,//2 have all moments 



s m {Hi) ■= \ x m dfii(x) 
Jo 

finite then so has /ii o // 2 and 

Sm (A*l ° /^) = S m (/ii) • S m (/i 2 ) 

for all m. 

If are absolutely continuous, with densities /1, / 2 respectively, then so is \i\ o /i 2 

and its density is given by the Mellin convolution: 



POD 

(/io/ 2 )(x):= / h{x/y)f 2 {y) [ 
Jo 



dy_ 

We will need the following modified beta distributions: 
Lemma 1.1. Letu,v,l>0. Then 

T(u + n/l)T(u + v) 



T(u + v + n/l)T(u) J n=0 
is the moment sequence of the probability measure 

(12) b(u + v,u,l) := — — — — -x lu ~^ (1 -x l Y~ l dx 

d(u, V) 

on [0, 1], where B is the Euler beta function. 
Proof. Using the substitution t = x l we obtain: 

T(u + n/l)T(u + v) B(u + n/l,v) 1 



T(u + v + n/l)T(u) B(u,v) B(u,v) J 

I 



[ e+^-^i-ty^dt 

Jo 



l — f 1 x^-'ii-xr-'dx. 



B(u, .... 

□ 

Note that if X is a positive random variable whose distribution has density f(x) and 
if / > then the distribution of X x l l has density lx l ^ 1 f(x l ). In particular, if X has beta 
distribution b(u + v, u, 1) then X x l l has distribution b(u + v, u, I). 

For u, I > we also define 

(13) b(u,u,l) := 



4 wojciech mlotkowski, karol a. penson, karol zyczkowski 

2. Applying Mellin convolution 

From now on we assume that p > 1 is a rational number, say p = k/l, with 1 < I < k, 
and that < r < p. We will show, that then the sequence A m (p, r) is a moment sequence 
of a probability measure fi(p,r), which can be represented as Mellin convolution of 
modified beta distributions. In particular, r) is absolutely continuous and we will 
denote its density by W p>r . The case when p is an integer was studied in [H 

First we need to express the numbers A m (p,r) in a special form. 

Lemma 2.1. If p = k/l, where k, I are integers, 1 < I < k and < r < p then 



(14) A m {p,r) = —==== ^- J ^(p)™ 



'2ixkl{k -I) \P-1J nj=i F («i + WO 

where c(p) = pP(p — l) 1_p , 

(15) dj- 

(16) & 

Proof. First we write: 

fmp + r\ r rT(mp + r) 




m J mp + r r(m + l)r(mp — m + r + 1) 
Now we apply the Gauss's multiplication formula: 

T{nz) = (2n)^n^V(z)T [z + ±) F (z + . . . F (z + ^ 



to get: 

r(m P + r) = r (*; (™ + £) ) = (2vr )(i^)/^^+-V2 Y[r(j + r + 3 ~ 1 
r(m + 1) = r = (27r)( 1 - i )/ 2 r +1 / 2 fj r (y + { 



and 



Y{mp — m + r + l) = T[{k — — + - — - 



k , 

(2nf- k+l ^ 2 (k - l^m(k-l)/l+r+l/2 J| r / ™ + 

j=l+l ^ 



m i r + j — I 
k-l 



It remains to apply them to (17). □ 
Now we need to modify enumeration of a's. 
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Lemma 2.2. For 1 < i < I + 1 denote 

(i-l)Jfe 



1. 



where [-J i/ie /Zoor function, so that 

l=ji<j 2 <---<Ji<k<k + l= 

For 1 < j < k define 



(18) 



I 

r+j 



k-l 



if 3 = ji, 1 < i < I, 
ifji <J <3i+i- 



Then the sequence {%} J=1 is a rearrangement of {ctj} ' =1 . Moreover, ifO < r < p = k/l 
then we have f3j < Sj /or a// j < fc. 

Proof. It is easy to verify the first statement. 

Assume that j = ji for some i < I. Then we have to show that 

r + ji-1 



which is equivalent to 



Ir + 1 



k 

k{i 



< 



I 



< ki. 



and the latter is a consequence of the fact that \_x\ < x and the assumption r < p = k/l. 
Now assume that ji < j < ji + \. We ought to show that 

r + j — 1 r + j — i 
k - k-l ' 

which is equivalent to 

Ir + Ij + k - / - ki > 0. 
Using the inequality \_x\ + 1 > x we obtain 

Ij + k — l — ki> l(ji + 1) + k - I - ki = Iji + k - ki > k(i - 1) + k - ki = 0, 
which completes the proof, as r > 0. □ 
Now we are ready to prove the main theorem of this section. 

Theorem 2.3. Suppose that p = k/l, where k,l are integers such that 1 < I < k, and 
that r is a real number, < r < p. Then there exists a unique probability measure 
fx(p, r) such that |ip is its moment sequence. Moreover /i(p, r) can be represented as the 
following Mellin convolution: 



fi(p, r) = b(«i, ft, I) o . . . o b(a k , 0k, ° 5 c (p), 



where 



c(p) := 



(p_ 

Consequently, ji(j>, r) zs absolutely continuous and its support is [0, c 



() 
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Note that the representation of densities in the form of Mellin convolution of modified 
beta distributions was used in different context in [7], see its Appendix A. 

Example. For the Marchenko-Pastur measure we get the following decomposition: 

(19) /i(2,l)=b(l,l/2,l)ob(2,l,l)o5 4 , 

where b(l, 1/2, 1) has density l/(iry/x — x 2 ) on [0, 1], the arcsine distribution with the 
moment sequence ( 2 f ^)4~ m , and b(2, 1,1) is the Lebesgue measure on [0,1] with the 
moment sequence l/(m+ 1). 



Proof. In view of Lemma 2.1 and Lemma 2.2 we can write 

k 



fJl r ("j + "V0 r (/3j) 



for some constant D. Taking m = we see that D — 1. □ 



Note that a part of the theorem illustrates a result of Kargin [12], who proved that 
if is a compactly supported probability measure on [0, oo), with expectation 1 and 
variance V, and if L n denotes the supremum of the support of the multiplicative free 
convolution power fi mn , then 

(20) lim — = eV, 

n— s>oo Tl 

where e = 2.71 ... is the Euler's number. The Marchenko-Pastur measure /i(2, 1) has 
expectation and variance equal to 1 and /i(2, l) Kn = fi(n + 1, 1), so in this case L n = 
(n + l) n+1 /n n (this was also proved in |24] and [TO]) and (20) holds. 

The density function for /i(p, r) will be denoted by W p>r (x). Since A m (p,p) = 
A m+1 (p, 1), we have 

(21) W p , p (x) = x-W Pil (x), 
for example 

(22) W 2 , 2 (x) = ^-y/x(4-x) on [0,4], 

which is the famous semicircle Wigner distribution with radius 2, centered at x = 2. 
Now we can reprove the main result of [13] . 

Theorem 2.4. Suppose that p, r are real numbers such that p > 1 and < r < p. Then 
there exists a unique probability measure /i(p, r), with support contained in [0, c(p)], such 
that {A m (p, r)}^ =0 is its moment sequence. 

Proof. It follows from the fact that the class of positive definite sequence is closed under 
pointwise limits. □ 



3. Applying Meijer G-function 

The aim of this section is to describe the density function W Ptr (x) of [i(p, r) in terms of 
the Meijer G-function (see [16] for example) and consequently, as a linear combination 
of generalized hypergeometric functions. We will see that in some particular cases W p ^ r 
can be represented as an elementary function. 
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Theorem 3.1. Let p = k/l > 1, where k, I are integers such that 1 < I < k, and let r be 
a positive real number, r < p. Then the density W P)T of the probability measure fi(p, r) 
can be expressed as 



(23) 



Wp, r {x) 



rp 



(p-l) r+1 / 2 V2k 



G 



k,0 / Qti, 



X 



71 



k.k 



Oik 



X 



c(p)> 



x G (0,c(p)), where c(p) =p p (p— l) 1 p and the parameters otj,(3j are given by (15) and 



Proof. Define 



rT(ap — p + r) 



rp,rV " ' T(a)T(ap-a-p + r + 2) 
If m is a natural number then 



px 



m + 1; 



mp + r 



r 



m J mp + r 



so (j)p^ is the Mellin transform of the density function W v<r of /i(p, r): 



In order to reconstruct W p r we apply the inverse Mellin transform: 



W P;r (x) 



1 

2vri 



C+ioo 



*> P)T {p) da, 



C—\oo 



PX 



see [31 HSl [19] for details. Putting m = a — 1 in ( 14 ) we get 

r( P - iy- l - r uU T Wj - l ' 1 + °^ 

pp- r y/2nkl(k-l) n- =1 rfo - 1// + a//) 

Writing the right hand side as $(a/l)c(p) tT , using the substitution a = lu and the 
definition of the Meijer G-function (see [16] for example) we obtain 



c { P y. 



W p>r (x) 



2ni 



C+ioo 



${<j/l)c{p) a x- a do- 



C—ioo 



I 

2~TXl 



C+ioo 
C—ioo 



(x l /c(p) 1 ) u du 



r(p-iy- r - 3/2 ck ,o 



a 1 ,...,a k 



j • • • j Hk 



Pi 



Z 



where z = x l /c(p) 1 , a- = a- — 1/7, /3- = /3j — 1/1. Finally we use formula (16.19.2) in 



[16J and obtain 
(24) 



r (p - l)P-r-3/2 



z 1 / l pP- r \ / 2~7ik 



■ ,0£k 



which is equivalent to (23). 



□ 



Now applying Slater's theorem (see (16.17.2) in [16]) we can represent W p>r as a linear 
combination of hypergeometric functions. 
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Theorem 3.2. For for p = k/l, with l<l<k,r>0 and x G (0,c(p)) we have 

(25) W p , r (x) ^t^^)^^^,^)^/^^^^) z\ z ^m-i/i 

h=l 

where z = x l /c(p) 1 , 

r (r) _ 1 \p-r-3/2 

(26) 7 (M,r)- U ' 



(27) c(h,k,l,r) - , ^ 

and ine parameter vectors of the hypergeometric functions are 

The most tractable case is p = 2. 

Corollary 3.3. For p = 2, < r < 2, the density function is 

sin ( r • arccos \/x/A ) 
(30) W 2 , r (x) = V nxl . r/2 

x G (0,4). In particular for r = 1/2 we have 



(3D »W*> = 

and /or r = 3/2 

(32) -WW ^ 

Note that if r > 2 then W 2 , r (x) < for some values of x G (0, 4). 

Proof. We take k = 2, I = 1 so that c(2) = 4, z = x/4 and 7(2, l,r) = r2 r /(8y/n). 
Using the Euler's reflection formula and the identity T(l + r/2) = T(r/2)r/2 we get 

T(l/2) 2sin(7rr/2) 



c(l,2,l,r) 
c(2,2,l,r) = 



r(l-r/2)r(l + r/2) 

r(-l/2) -2cos(7rr/2) 



r((l-r)/2)r((l + r)/2) A 
We also need formulas for two hypergeometric functions, namely 

r — r 1 
2' ~2~' 2 
1 + r 1 -r 3 



2-^1 

/I 

2^1 



2 ' 2 ' 2 



cos(r arcsin \/z), 
sin(r arcsin y/z) 



Tyfz 
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see 15.4.12 and 15.4.16 in [16]. Now we can write 

sin(7rr/2) cos ( r arcsin a/x/4) — cos(7rr/2) sin (r arcsin y/x/A 
W-jA-r) V 



7TX 



l-r/2 



sin {^ir/2 — r arcsin y/x/Aj sin arccos y/x/A 



7TX 1 r l 2 TlX 1 r l 2 

which concludes the proof. □ 
Remark. Note that 

W2.1 (Vg 1 (X 

(33) = fTw ^ 

It means that if X,Y are random variables such that X ~ p(2, 1) and V ~ /x(2, 1/2) 
then X 2 ~ AY. This can be also derived from the relation A m {2, l/2)4 m = A 2m (2, 1). 

4. Some particular cases 

In this part we will see that for k = 3 some densities still can be represented as 
elementary functions. We will need two families of formulas (cf. 15.4.17 in |16j). 

Lemma 4.1. For c ^ 0, —1, —2, . . . we have 

(34) 2 F 1 ^,^; C L")=2 c - 1 (l + Vr^) 1 - c , 



(35) 2 Fi 



fc+1 c-2 \ 2 C_1 / , u-c/ / x 

( — ,— ; c = -—(l + vT^) 1 c (c-l + v / T^). 

Proof. We know that 2 .Fi(a, 6; c| z) is the unique function / which is analytic at z = 0, 
with /(0) = 1, and satisfies the hypergeometric equation: 

z(l - z)f"{z) + [c - (a + b + l)z] f'(z) - abf(z) = 

(see (3]). Now one can check that this equation is satisfied by the right hand sides of 
(34) and (j35l for given parameters a, b, c. □ 



Now consider p = 3/2. 
Theorem 4.2. Assume that p = 3/2. Then for r = 1/2 we have 

/ , \2/3 / n2/3 

(1 + v/1 -4x 2 /27) - (1 - y/l -4x 2 /27) 

(36) WW*) = " 25/33-1/2^2/3 — ' 

/or r = 1 : 

,37) W 3/2 , 1{X) = + - (> - " 



2 4 /%x 1 /3 

/ , \2/3 / n2/3 

1/2 1/3 I 1 + V 71 " 4 *727j - (l - a/1 - 4x 2 /27j 



+3 1 /V/ 3 - 

2 5 / 3 7T 

and, finally, Wy^^x) = x ■ Ws/2,i(x), with x G (0,3-^/3/2). 
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Proof. For arbitrary r we have 

2 1 - 2r / 3 sin (2vrr/3) / 3 + 2r r -2r 2 1 

W 3 /2,r{x) - Q3 /2-r„ * F * [ ~ ' 3' 3' 3 



33/2-r^ 
2 (4-2r)/3 r sin ((l _2r)7T/3) 



r/3-1/2 



33/2-r^ 



3-^2 



5 + 2r 1 + r 1 - 2r 4 2 



6 ' 3 ' 3 ' 3' 3 



z z 



(r+l)/3-l/2 



r(l + 2r) sin ((1 + 2r)vr/3) / 7 + 2r 2 + r 2-2r 5 4 

— 3-^2 



2 (l+2r)/3 3 3/2-r 7r 



6 ' 3 ' 3 



3' 3 



z z 



(r+2)/3-l/2 



where z = 4x 2 /27. Ifr = l/2orr = l then one term vanishes and in the two others 
the hypergeometric functions reduce to 2-^1- 



For r = 1/2 we apply (34) to obtain: 



W 3 /2,i/2(x) = 

z -l/3 



2 1 /331/2 7T 



g-^ 3 /l -1 1 

2 1 /3 3 i/2 7r 2 1 U'T' 3 



2- 2 /3 (l + v ^3^) 2 / 3 _ 



4/3 



2 5/3 3 l/2 



7T 



2^1 



5 15 
6'3' 3 



,1/3 



2 5 /33l/2 



7T 



2 2/3 (! + v ^3^) 



-2/3 



-1/3 



2 • 3V27T 

Z^/ 3 
= 2 ■ 3*/% 



(i + vtt^) 



2/3 2 1/3 /l-^T 3 



(l + v/T 3 ^) 



2/3 Z 



2 ■ 3*/ 2 7T 
1/3 



2 • 3V2tt 



2/3 



2/3 



and this yields (36). 



For r = 1 we use (35): 



z ~ x l % /5 — 2 2 \ 

w y^ = vJH 2F \6>T ; 3 z ) 



z 1 / 6 (7 -I 4 
+ 2V% 2 H 6'X ! 3 



-1/6 



(i + v^^fc/T^-lJ + ^l + vT 



,1/6 



,-1/6 



(i + v / T^7) 1/3 (3vT^i-i) 



47T 

Now we have 



-1/6 



-1/3 



(i_ v /T37 ) 1 / 3 ( 3> /r=I+i) . 



(3\/r^ + 1) 



47T 



(1 + v/r^) 173 (sv/r^ - 1) = - (1 + Vi^y* (3 - sv/r^ - 2) 

= -3Z 1 / 3 (1 - VT^~zf 3 + 2 (1 + VT^) 1/3 

and similarly 

(1 _ v/T37) 1/3 (3v / T^+ 1) = 3^ 3 (1 + VT^~zf 3 - 2 (1 - s/l^z) 
Therefore 



1/3 



1/3 



W 3 / 2 ,l {x) 



,-1/6 



(i + v / r^) 1/3 -(i-v / r^) 



1/3 



+ 



3^/6 



47T 



((i + v / r^) 2/3 -(i-v / r^) 



2/3 



which entails (37). The last statement is a consequence of (21). 



□ 
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The dilation D2/i(3/2, 1/2), with the density Ws/2,i/2{x/2)/2, is known as the Bures 
distribution, see (4.4) in [21J . Moreover, the integer sequence A m A m (3/2, 1/2), moments 
of D 4 /i(3/2, 1/2), appears as A078531 in [20] and counts the number of symmetric 
noncrossing connected graphs on 2n + 1 equidistant nodes on a circle. The axis of 
symmetry is a diameter of a circle passing through a given node, see [6]. 

The measure /i(3/2, 1) is equal to /i(2, l) m / 2 ; the multiplicative free square root of 
the Marchenko-Pastur distribution. 

For the sake of completeness we also include the cases p — 3, r — 1 and p = 3, r = 2, 
which have already appeared in [XTl | 



Theorem 4.3. Assume that p = 3. T/ien /or r = 1 we have 

2/3 



(38) 

for r = 2: 
(39) 



W 3 ,i(x) 



3 (l + v 7 ! -4x/27)' 



22/3^1/3 



24/331/2^2/3 ^ + ^1-4^/27) 



1/3' 



W 3 ,2(X) 



9 (l + Ax/27 



4/3 



24/3^2/3 



_\ 2/3 



25/333/2^1/3 ^ + ^1-4^/27) 

and, finally, W 3j3 (x) = x ■ ^3,1(2;), with x G (0,27/4). 
Proof. For arbitrary r we have 

2( 6 ~ 2r )/ 3 sin (vrr/3) /r 3-r -r 2 I 

W^3,r(x) = 3^2 

2 (4-2r)/3 r sin ((i _2r)7r/3) 

3 3 - r 7T 

r(r — 1) sin ((1 — r)n/3) 
2( 1+2r )/ 3 3 3 - r 7r 



3^2 



3-^2 



3' 6 ' 6 ' 3' 3 

1 + r 5 - r 2 - r 4 2 
3 ' 6 ' 6 ' 3' 3 

2 + r 7-r 4-r 5 4 



z \z 



(r-3)/3 



z 2; 



(r-2)/3 



6 



6 



' 3'3 



2 z 



(r-l)/3 



where z = Ax/21. For r = 1 and r = 2 we have similar reduction as in the previous 
proof. Here we will be using only (34). 
Taking r = 1 we get 

2 l/3 z -2/3 



W31 X 



3 3 / 2 7T 
,-2/3 



1 -1 2 

3'~6~' 3 



,-1/3 



2 1 /333/2 ?r 



iF 1 



33/2 



71 



(1 + v / r^) 



1/3 z 



-1/3 



3 3 / 2 7T 
2 / 3 -.1/3 



2 1 4 

3' 6' 3 



-1/3 



33/2^^2/3 ( 1 + v ^^) 1 /3' 



which implies (38). 
Now we take r = 2: 
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(l + v/r^) 473 -^ 3 



2 ■ ^/Hz 1 ^ (1 + VT^~z) 



2/3 : 



and this gives us (39). Finally we apply rt2JJ). □ 



Note that the measure /i(3, 1) is equal to /x(2, 1) K2 , the multiplicative free square of 
the Marchenko-Pastur distribution. 

5. Graphical representation of selected cases 



The explicit form of W p%r (x) given in Theorem 3.2 permits a graphical visualization for 
any rational p > and arbitrary r > 0. We shall represent some selected cases in Figures 
1-9. These graphs which are partly negative are drawn as dashed curves. In Fig. [I] the 
graphs of the functions Wz/2,r{%) for values of r ranging from 0.9 to 2.3 are given. For 
r < 3/2 these functions are positive, otherwise they develop a negative part. In Fig. 2 
we represent W§/2, r ( x ) for r ranging from 2 to 3.4. In Fig. 3 we display the densities 
W PjP (x) for p = 6/5,5/4,4/3 and 3/2. All these densities are unimodal and vanish at 
the extremities of their supports. They can be therefore considered as generalizations 



of the Wigner's semicircle distribution ^2,2(^)5 see equation (22). In Fig. 4 we depict 
the functions W / 4/3, r (^), for values r ranging from 0.8 to 2.4. Here for r > 1.4 negative 
contributions clearly appear. In Fig. 5 and 6 we present six densities expressible through 



elementary functions, namely W^/2,r{ x ) for r = 1/2, 1, 3/2, see Theorem 4.2 and W^^ r {x) 



for r = 1, 2, 3, see Theorem 4.3, In Fig. 7 the set of densities W Pt i(x) for five fractional 
values of p is presented. The appearance of maximum near x — 1 corresponds to the 
fact that p,(p, 1) weakly converges to Si as p — > 1 + . In Fig. 8 the fine details of densities 
W p> i(x) for p = 5/2, 7/3, 9/4, 11/5, on a narrower range 2 < x < 4.5 are presented. In 
Fig. 9 we display the densities W P) i(x) for p = 2,5/2,3,7/2,4, near the upper edge of 
their respective supports, for 3 < x < 9.5. 

The Fig. 10 summarizes our results in the p > 0,r > quadrant of the (p,r) plane, 
describing the Raney numbers (c.f. Fig. 5.1 in [13j and Fig. 7 in [18J). The shaded 
region E indicates the probability measures ^(p,r) (i.e. where W Pir (x) is a nonegative 
function). The vertical line p = 2 and the stars indicate the pairs (p, r) for which 
W P7r (x) is an elementary function, see Corollary 3.3, Theorem 4.2 and Theorem 4.3 



The points (3/2, 1) and (3, 1) correspond to the multiplicative free powers MP m / 2 and 
MP m of the Marchenko-Pastur distribution MP. Symbol B at (3/2, 1/2) indicates the 
Bures distribution and SC at (2, 2) denotes the semicircle law centered at x = 2, with 
radius 2. 

It is our pleasure to thank M. Bozejko, Z. Burda, K. Gorska, I. Nechita and M. A. Nowak 
for fruitful interactions. 
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Figure 1. Raney distributions W^/2, r { x ) with values of the parameter r 
labeling each curve. For r > p solutions drawn with dashed lines are not 
positive. 
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X 

Figure 2. As in Fig. [I] for Raney distributions W 5 /z >r (x). 
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Figure 5. Raney distributions Wz/2,ri x ) with values of the parameter r 
labeling each curve. The case Wz/2,i{ x ) represents the multiplicative free 
square root of the Marchenko Pastur distribution. 
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Figure 6. Raney distributions W% yr {x) with values of the parameter r 
labeling each curve. The case (x) represents the multiplicative free 
square of the Marchenko Pastur distribution. 
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Figure 7. Raney distributions W Pt i(x) with values of the parameter p 
labeling each curve. The case Wz/2,i{ x ) represents the multiplicative free 
square root of the Marchenko-Pastur distribution, MP m / 2 . 
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Figure 8. Tails of the Raney distributions W Pt i(x) with values of the 
parameter p labeling each curve. 
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Figure 9. As in Fig. [8] for larger values of the parameter p. 




X 



Figure 10. Parameter plane (p,r) describing the Raney numbers. The 
shaded set S corresponds to nonnegative probability measures r). 
The vertical line p = 2 and the stars represent values of parameters 
for which W p>r (x) is an elementary function. Here MP denotes the 
Marchenko-Pastur distribution, MP Ks its s-th free mutiplicative power, 
B-the Bures distribution while SC denotes the semicircle law. For p > 1 
the points (p,p) on the upper edge of £ represent the generalizations of 
the Wigner semicircle law, see Fig. [3} 
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